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Abstract 

The interplay of the effects of geometry and collective motion on d-a correlation functions is 
investigated for central Xe+Au collisions at E/A=50 MeV. The data cannot be explained without 
collective motion, which could be partly along the beam axis. A semi-quantitative description of 
the data can be obtained using a Monte-Carlo model, where thermal emission is superimposed on 
collective motion. Both the emission volume and the competition between the thermal and collec- 
tive motion influence significantly the shape of the correlation function, motivating new strategies 
for extending intensity interferometry studies to massive particles. 
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Nuclear collisions provide the only means to study highly excited nuclear matter and 
its phase transitions under laboratory-controlled conditions. Such collisions produce highly 
excited nuclear systems that persist momentarily. Space-time information about the fate of 
these systems can be accessed through the dependence of two-particle correlation functions at 
low relative momentum on Boson or Fermion symmetries and on the particle mutual nuclear 
and Coulomb interactions [1-4]. Particularly successful investigations using proton-proton 
correlations at intermediate energies and pion-pion correlations at high energies have been 
performed in the last decades [1-8] . Extending these studies to all particle species produced 
during a reaction represents an important objective because different particles may originate 
at different stages of the reaction where the densities are different. Correlation functions 
between complex particles and heavy fragments may be more relevant than light particle 
correlation functions to stages of the reaction where multifragmentation and the liquid-gas 
phase transition are expected to occur [8-14]. 

Since thermal velocities decrease with particle mass, v th oc (m) -1 ^ 2 , while collective ve- 
locities do not depend on mass, the latter become increasingly important as more massive 
particles are considered [15,16]. This influences the connection between the space-time struc- 
tures of the system and the measured correlations [15-17] for such particles. In this article, 
we explore the effects of collective motion on correlation functions constructed with pairs of 
complex particles interacting by means of the mutual nuclear and Coulomb final state inter- 
actions. We present our ideas by studying deuteron-alpha correlation functions measured 
in 129 Xe+ 197 Au collisions at E/A=50 MeV. We demonstrate that the correlation data can 
be understood only if one incorporates significant effects induced by collective motion and 
temperature along with the space-time geometry. 

We studied central 129 Xe+ 197 Au collisions at E/A=50 MeV, using Au targets of 3 mg/cm 2 
areal density and 129 Xe beams, having an intensity of about 108 pps, from the K1200 cy- 
clotron of the National Superconducting Cyclotron Laboratory at Michigan State University. 
Isotopically resolved particles with 1<Z<10 were detected with nine telescopes of the Large 
Area Silicon Strip Array (LASSA) [18,19]. Each telescope consisted of one 65 /im single- 
sided silicon strip detector, one 500 /im double-sided silicon strip detector and four 6-cm 
thick CsI(Tl) scintillators. The strips of the silicon detectors provided an angular resolution 
of about ±0.43°. The center of the LASSA device was located at a polar angle of 6=35° 
with respect to the beam axis, covering polar angles of 12° < 9 < 62° and azimuthal angles 
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24° < < 156°. Impact parameters were selected by the multiplicity of charged particles, 
measured with LASSA and the 188 plastic scintillator- CsI(Tl) phoswich detectors of the 
Miniball/Miniwall array [20]; the combined apparatus covered 80% of the total solid an- 
gle. Quasi-central collisions, with reduced impact parameters of the order of b/b max <0.3, 
were selected by a gate on the measured charged-particle multiplicity distribution. This re- 
duced impact parameter gate can contain non-negligible contaminations from mid-peripheral 
events. Statistics requirements characterizing two-particle correlation analyses prevented us 
from using more selective impact parameter filters as those described in ref. [21]. 

The data points on Fig. 1 show the d-a correlation function, 1 + R(q), defined in terms 
of the two-particle differential multiplicity, dM da (p d ,p a ) /dp d dp d , and the single particle 
differential multiplicities, dM d (p d ) jdp\ and dM a (p a ) jdp\: 



HdM da {p d , p a ) /dpldpl = C-[1 + R (q)) ■ Y,dM d (p d ) /dpi ■ ZdM a (p Q ) I dpi (1) 



where p d and p a are the laboratory momenta of the two coincident particles, q = /i • v re i 
is the momentum of relative motion, and C is a normalization constant chosen such that 
< R(q) >=0 for large g-values where final-state interaction effects are negligible. The sums 
on each side of Eq. ([I]) extend over all particle energies and angles contributing to each bin of 
q. The product of the single particle yields, on the r.h.s. of Eq. ([!]) has been constructed by 
mixing particles from different events [22]. We apply the same impact parameter selection 
to both single and two particle multiplicity distributions in Eq. ([I]). 

The d-a correlation function exhibits a large minimum at small relative momentum due 
to the mutual Coulomb repulsion. The structures at larger relative momentum arise from 
the d-a nuclear interaction. In particular, the sharp peak at q ^42 MeV/c corresponds to 
the first excited state of 6 Li at E*( 6 Li) = 2.186 MeV (J w =3+, T=0.024 MeV, r c /r=1.0) and 
the broad peak around 84 MeV/c stems mainly from the resonance at E*( 6 Li) = 4.31 MeV 
(J 7r =2 + , T=1.3 MeV, r c /r=0.97) with small contributions from the resonance at E*( 6 Li) 
= 5.65 MeV (J 7r =l + , T=1.9 MeV, r c /r=0.74). Theoretically, the correlation function is 
commonly calculated using the angle-averaged Koonin-Pratt equation [1,5-8]: 



The angle-averaged kernel K(q, r) is calculated from the radial part of the d-a elastic scat- 
tering relative wave function. The source function, S(r,q,P), is defined as the probability 
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for two-correlated particles with relative momentum q and total momentum P to be sep- 
arated by a distance r when the last of the two particles is emitted. Most studies using 
Eq. ([2]) have assumed that the source function depends only on the relative distance, r, 
and not on relative momentum, q [1,2,5-8,23,24]. Under this assumption, Gaussian source 
profiles, S(r) oc exp (— r 2 //^), depending only on the geometry of the system, have failed in 
reproducing the line-shape of the correlation functions [2,23,24]. In particular, theoretical 
correlation functions that reproduce the magnitude of the E*=2.186 MeV peak, typically 
over-predict the E* «4.31 MeV peak [24]. Even if one accounts for long-lived secondary de- 
cays of unstable fragments, the situation does not improve because the relative magnitude 
of the two peaks does not change. Indeed, as it was shown in the case of proton-proton 
correlation functions, long-lived secondary decays reduce the magnitude of the correlation 
peaks without changing their shape [7] . 

Stimulated by the difficulties encountered by previous approaches, we explore in details 
how the line-shape of a d-a correlation function can be affected by the interplay of the 
geometry of the source with thermal and collective motion. We introduce a Monte Carlo 
simulation where deuterons and alpha particles are emitted at random positions Td and f a 
within a spherical single-particle source having radius Rq. This source evaporates particles 
while undergoing collective expansion. Then the velocities of the two emitted particles con- 
sist of a thermal and a collective component, Vd = Vdjh+Vd,coii ifd) and v a = v at th+Va,coU (r a ). 
The thermal velocities, Vd,th an d v a jh, are randomly assigned according to Boltzmann distri- 
butions corresponding to an assigned temperature, T=4 MeV, consistent with measurements 
of isotopic and excited states thermometers [25] . The collective velocity components are pa- 
rameterized as vector fields, Vd,coii ifd) and v a ,coii depending on particle positions, 
and f a . Simulated particle trajectories are filtered through the acceptance of the LASSA 
array. Finally, the relative distance, r = fd — f a , and the momentum of relative motion, 
q '= fi(vd — v a ), are evaluated and used to construct the two-particle emitting source, S(r, q), 
by integrating over total momentum, P = \pd + p a \- The source function is then used in 
Eq. ^ to calculate the correlation function, R(q). 

We first compare our data to calculations performed with a pure-thermal emitting source 
where particle velocities have only their thermal components (corresponding to the assigned 
temperature, T=4 MeV) and no collective components (vd lCO ii — Va,coii = fy- The main impact 
of varying the radius of the source, R , is to change the widths of the peaks [7]. The width of 
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FIG. 1: d-a correlation function. Top panel: overall view of the correlation function. Bottom 
panel: view emphasizing the second resonance at q ~84 MeV/c. Circles: d-a correlation function 
measured in Xe+Au collisions at E/A=50 MeV. Lines: calculated correlation functions for a pure- 
thermal spherical source with i?o=10.5 fm, T=4 MeV (dotted line), and an ellipsoidal source with 
(i?z=16fm,i?x=8fm,i?y=8fm), T=4 MeV and collective motion parameter (3 co u=0.1 (solid line), 
0.2 (dashed line) and 0.3 (dot-dashed line). 
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the first sharp peak at q=42 MeV/c is dominated by the finite resolution of the LASSA array, 
leaving no sensitivity to the source size. In contrast, the width of the broad peak at 84 MeV/c 
is found to scale almost linearly with the radius of the source, i.e. Aq FWHM (MeV /c)~25- 
1.2-i?o(f m )- This sensitivity provides the opportunity to determine the source size from 
measured d-a correlation functions. Calculations with a spherical source with i?o=10.5 fm 
approximately reproduce the width of the peak at q ~84 MeV/c. As for the p-p correlation 
function, the main effect of long-lived secondary decays is to reduce the magnitude of the 
correlation peaks without changing their shape [7]. We take this into account by fitting the 
measured correlation function R exp (q) by Rf as t(<l) = A-.R(g), where A represents the fraction 
of d-A pairs from the fast source and R(q) is calculated from Eq. The dotted line on 
Fig. 1 shows the calculated correlation function obtained by choosing A=0.8. Choosing this 
value of the A parameter normalizes the correlation function to the first peak at 42 MeV/c, 
but over-predicts the height of the second peak at q ^84 MeV/c and under-predicts the 
valley between the two peaks, as seen in the lower panel of Fig. 1 showing an expanded 
view of R(q) for 50< q <130. These difficulties, already encountered in previous works and 
confirmed by our high resolution data, cannot be overcome by changing Ro, A, the source 
temperature or by altering the source geometry. Indeed, the d-a kernel function, K(r,q), 
has the same spatial extension close to the resonances at q ~42 and 84 MeV/c: any change 
in the geometry of the source would change the height of the two peaks by the same factor 
and leave their ratio unmodified. In the following, we focus our attention upon this long- 
standing and well documented [24] problem and defer optimizations with respect to source 
size, temperature and geometry to later investigations. 

In order to explain the observed discrepancies, we recognize that our selection of events 
with high measured charged-particle multiplicities with b < 0.3 includes events with non- 
negligible collective motion, both along and perpendicular to the beam axis [21,26,27]. For 
simplicity we parameterize collective velocities as Vd,coii ifd) = cj3 co ur d /R and v a)CO ii {r a ) = 
c^coiifa/ Ro , where c is the velocity of light and (3 co u represents the maximum collective 
velocity of a particle emitted at the surface of the system, i.e. when = r a = Ro. This 
parameterization recalls self-similar radial expansion profiles already used in the literature 
[21,26-29]. 

Choosing deuteron-alpha particle pairs of from the simulation that are detected in 
LASSA, we construct the source function S(r) that describes their relative separation at 
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emission and found that collective motion induces an apparent reduction of the source size, 
consistent with studies at higher incident energies [16,17]. The thick solid line in Fig. 2 
shows the profile of S(r) for a source with (3 co u=0.2, i? =10.5 fm and T=4 MeV; the dotted 
line shows S(r) for a pure thermal source with /3 co ii=0, Rq=10.5 fm and T=4 MeV. In this 
pure-thermal case, particles have velocities that are independent of the position where they 
are emitted. Since particles emitted from any location in the source have an equal probabil- 
ity of being detected by the LASSA array, the corresponding correlation function probes the 
whole spatial extent of the source. If collective motion is present, however, particles emitted 
from regions where the collective motion points away from LASSA will be included in the 
correlation function only if their thermal velocity is larger than the local collective velocity. 
This limits the portion of the emitting source that can be probed by the experimental data 
and the source size appears reduced. In the particular case shown on Fig. 2, the apparent 
source size is reduced by as much as 50%. 

The source size reduction is just one aspect of the influence of collective motion on 
correlations. Another is a dependence of the correlation function on the total momentum, 
P = pi + f>2 of the particle pairs [1,2-8,16,17]. Such total momentum-dependence has also 
been interpreted at intermediate energies as a consequence of emission from cooling and 
expanding systems [1-4,23]. In Fig. 2 we show that the extracted source size for particles 
within a small relative momentum gate, g=0-50 MeV/c, is even smaller (thick dashed line) 
than the size of source corresponding to all extracted d-a pairs (solid line). Thus, collective 
motion makes the source function both r- and g-dependent. It happens because the d- 
a relative momentum, q, has thermal and collective components, q = cfth + Qcoii (t)- The 
collective component, given by q co u (f) = // C| ^°" r , depends on the relative distance, making 
the source both r- and g-dependent, i.e. S = S(r,q), with small relative momentum values 
being positively correlated with correspondingly small relative distances. This correlation 
between relative distance, r, and relative momentum, q, is an unexplored effect of collective 
motion on two-particle emitting sources. In the case of a pure-thermal model (/3 C o«=0), the 
relative momentum is independent of r and contains only a thermal component, q ~ q t h = 
A* (v<i,th ~ Va,th)- Indeed, in Fig. 2 we show that the pure-thermal source profile obtained by 
gating at small relative momentum, g=0-50 MeV/c (thin dashed line) is exactly the same 
as the original ungated one (dotted line). 

To show how collective motion influences the shape of the d-a correlation function, we 
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FIG. 2: Calculated d-a sources for a spherical single-particle source with Ro=10.5 fm and T=4 
MeV. Pure-thermal source calculations, f3 co u=0.0, are represented by a dotted line (all qs) and by 
a thin dashed line ((7=0-50 MeV/c). The thick solid and thick dashed lines represent the spherical 
source with /3 CO //=0.2 and, respectively, without and with a g=0-50 MeV/c gate applied. 

separate it into two components, R(q) = Rc(q) + RN(q), where Rc{q) and Rn(<i) correspond 
to the correlation functions generated separately by the mutual d-a Coulomb and nuclear 
interactions, respectively. The left and the right panels of Fig. 3 show, respectively, 1+R c (q) 
and l + R N (q) obtained for a source with _R =10.5 fm, T=4 MeV, and (3 co u=0.0 (dotted line), 
f3 co u=0.l (solid line), /3 co u=0.2 (dashed line) and /3 co u=0.3 (dot-dashed line). The nuclear 
correlation functions (right panel) are all normalized to the integral of the first peak at g=42 
MeV/c, not shown on Fig. 3. 

The Coulomb correlation function (left panel, Fig. 3) displays one prominent anti- 
correlation at small relative momentum due to the mutual Coulomb repulsion between the 
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FIG. 3: Effects of collective motion on Coulomb (left panel) and Nuclear (right panel) correla- 
tion functions. Dotted, solid, dashed and dot-dashed lines correspond to calculations performed, 
respectively, with /3 coll =0.0, /3 coa =0.1, /3 co u=0.2 and (3 co u=0.3. 

deuteron and the alpha particle. The width of the Coulomb dip increases with collective 
velocity, consistent with a decreasing apparent source size caused by the collective motion in 
the system. This result shows that the possible existence of collective motion must be taken 
into account in order to extract source sizes and emission times from Coulomb dominated 
correlation functions [30]. 

On the right panel of Fig. 3 we show the calculated nuclear correlation functions in the 
region where the second peak is observed (q— 50-110 MeV/c). The magnitude of the peak 
at 84 MeV/c decreases with increasing (3 co u. In order to understand this effect, we point out 
that the relative momentum regions where resonant peaks exist, are mostly dominated by 
d-a pairs emitted at very small relative distances, r <2 fm, corresponding to the average 
d-a spatial separation in the resonance. Under these conditions, the two coincident parti- 
cles have nearly the same collective velocities because they originate from nearly the same 
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location, resulting in a collective component of relative momentum given by q co u (r rs 0) ~ 
and q ~ g t /j. Then, the relative momentum distribution can be described as a Boltzman 
distribution at the local temperature, T, of the source. The uncorrelated two-particle spec- 
trum in the denominator of Eq. ([!]), on the other hand, is calculated by selecting deuterons 
and alpha particles detected in different events. Either particle can originate from anywhere 
in the source, thus the collective velocities of d and a can be very different and q co u ^ 0. 
Then the uncorrelated relative momentum spectrum is nearly exponential but has a larger 
slope, T MIX > T, as compared to that of the coincidence two-particle spectrum. The re- 
sulting correlation function, obtained by folding the coincidence spectrum with the kernel 
K (r, q) and by dividing by the uncorrelated spectrum, displays an exponential fall-off in 
those relative momentum regions where resonances are observed. The slope, T app , of such 
fall-off is approximately given by the relation 1/T app — 1/T — 1/Tmix and it is the origin of 
the attenuation of the height of the second peak observed in the right panel of Fig. 3. This 
effect becomes increasingly important as the masses of the particles are increased, reflecting 
the more relevant role played by collective motion in the case of heavier particles. Previous 
studies with complex particles have commonly assumed that the temperature dependence 
in the numerator and denominator of Eq. Q are the same and cancel out in calculating 
their ratio [23]. As a consequence, the resulting correlation function resulted sensitive only 
to the volume of the source and not to the temperature [23]. In contrast, we show on Fig. 
3 that even small amounts of collective motion make the effective temperature, Tmix, of 
the denominator larger that the thermal temperature. The resulting correlation function is 
then sensitive not only to the volume of the source but also to its temperature through the 
interplay of geometry with collective motion. Our result provides the basis for unified ap- 
proaches to correlation function analyses to extract both space-time and thermal properties 
of emitting sources. 

To compare calculated correlation functions to our data in Fig. 1, we make more realistic 
assumptions about the geometry and the collective motion of the emitting source. Trans- 
verse collective expansions of about 2 MeV per nucleon have been previously observed in 
central Xe+Sn symmetric collisions at E/A=50 MeV by studying complex fragment kinetic 
energy spectra [21]. It has also been observed that for central events selected by high mea- 
sured charged particle multiplicities, the longitudinal collective motion in more asymmetric 
Xe+Au reactions is larger [31]. Thus, we simulate the emission of deuterons and alpha 
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particles to occur from an ellipsoidal source characterized by radii Rz = 2Rx = 2Ry, where 
Z indicates the beam direction. Dcutcrons and a particle collective velocities are described 
as v d>co u (f d ) = cficou (f d /R z ) and v^ coU (f a ) = c(3 coU (f a /R z ), respectively. This parameteri- 
zation simulates the expected radial collective motion in the center of mass of the reaction 
and provides collective velocity components to those particles emitted close to the edges 
ij'd^'a ~ Rz) of the ellipsoid, representing the effects of incomplete stopping of the incident 
momenta. 

Correlations functions have been calculated for a ranges of radii, (Rz, Rx, Ry), and for 
different values of the radial velocity parameter, (3 co u=0.l, 0.2 and 0.3. We searched for those 
combinations of (Rz,Rx,Ry) and f3 co u that provided a width of the peak at 84 MeV/c in 
reasonable agreement with the data. A reasonable description of that width is obtained for 
R z =l5-20 fm, R x =7.5-10 fm, i2 y =7.5-10 fm and (3 coU wO.1-0.2. In the bottom panel of 
Fig. 1, we only show the results corresponding to R Z =IQ fm, R x =Ry=8 fm. The different 
curves correspond to (3 co u=0.1 (solid line), (3 co u=0.2 (dashed line) and (3 co u=0.3 (dot-dashed 
line). All simulated correlation functions shown in Fig. 1 are normalized to reproduce 
the magnitude of the first peak at g=42 MeV/c. Values of /3 CO «=0. 1-0.2, providing a bt- 
ter agreement with the data, induce expansion velocities in the central region of the order 
of Prad ~ Pcoii/2 ~ 0.05-0.1, consistent with previous measurements of transverse flow for 
Xe+Sn and Xe+Au collisions [21,31]. Maximum collective velocities at the edges of the 
source (r = R z ) are of the order of (3 co u ?s0. 1-0.2, consistent with the velocity of projectile- 
like and target-like remnants in the center of mass of the studied reactions. By introducing 
collective motion (solid line), we find that one can reduce significantly the observed discrep- 
ancies between experimental data and the pure-thermal source calculation (dotted line). 
The height of the peak at 84 MeV/c still remains slightly over-predicted. The valley be- 
tween the peak at 42 and 84 MeV/c is well described by the calculations. More refined 
descriptions of collective motion (expansion, rotation, etc.) and of the break-up geometry 
might provide somewhat better agreement. However, our main goal is to demonstrate that 
the interplay between thermal and collective motion introduces a g-dependence in the source 
function that can account for the long-standing difficulties in describing the line-shape of 
d-a correlation functions. Given the range of radii that yield reasonable descriptions of the 
data and assuming that the source contains about 80% of the total mass [21,31], we find that 
the source volume is consistent with average nuclear densities ranging between p ^0.2p and 
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0.4p (with p =0.16 fm~ 3 , saturation density). Such range is consistent with the freeze-out 
density values commonly employed in statistical models of multifragmentation [32] . 

The present work is limited to the study of inclusive d-a pairs, without any condition 
imposed on their velocities or total momentum. Such inclusive data contain contributions 
from different stages of the reaction. The existence of these multiple sources of deuteron and 
a particles has been predicted by the Expanding Emitting Source model [33] and is expected 
to influence the shape of correlation functions [7]. We expect that a quantitative comparison 
to experimental data can be improved by implementing our simplified single-source model 
within a more realistic scenario where multiple sources with different timescales are included. 
Such ideas can be also tested by gating on the total momentum of the correlated particles. 

In summary, in this work we show that the interplay between source geometry and collec- 
tive motion plays a key role in determining the line-shape of correlation functions constructed 
with complex particles. In addition to the previously studied effects induced by collective 
motion, we observe also the presence of a significant relative momentum dependence in the 
emitting source. We apply our ideas to study a d-a correlation function measured in Xe+Au 
collisions at E/A=50 MeV. Realistic assumptions about the volume, the temperature and 
the collective motion in the system allow us to obtain a semi-quantitative description of the 
data and explain the difficulties encountered in previous studies. The presented results also 
show, for the first time, that the detailed line-shape of high resolution correlation functions 
probe not only the space-time geometry of the emitting source but also its interplay with 
thermal and collective motion. 
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